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ABSTRACT 

The model of stochastic acceleration of particles by turbulence has been successful in explaining 
many observed features of solar flares. Here we demonstrate a new method to obtain the accelerated 
electron spectrum and important acceleration model parameters from the high resolution hard X-ray 
observations provided by the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI). In 
our model, electrons accelerated at or very near the loop top produce thin target bremsstrahlung 
emission there and then escape downward producing thick target emission at the loop footpoints. 
Based on the electron flux spectral images obtained by the regularized inversion of the RHESSI 
count visibilities, we derive several important parameters for the acceleration model. We apply this 
procedure to the 2003 November 03 solar flare, which shows a loop top source up to 100-150 keV in 
hard X-ray with a relatively flat spectrum in addition to two footpoint sources. The results imply 
presence of strong scattering and a high density of turbulence energy with a steep spectrum in the 
acceleration region. 

Subject headings: acceleration of particles — Sun: flares — Sun: X-rays, gamma rays 



1. INTRODUCTION 

It is well established that the impulsive phase hard 
X-ray (HXR) emission of solar flares is produced by 
bremsstrahlung of nonthermal electrons spiraling down 
the flare loop whil e losing energy primarily via elastic 
Coul omb collisions (|Brownlll97lHHud son 1972; Pc trosianl 
Il973h . Thus, HXR observations provide the most direct 
information on the spectrum of the radiating electrons 
and perhaps on the mechanism responsible for their ac- 
celeration. The common practice to extract this infor- 
mation has been to use the parametric forward fitting 
of HXR spectra to emission by an assumed spectrum, 
usually a power-law with breaks and cutoffs (or plus 
a thermal co mponent), of the rad iating or accelerated 
electrons fe.g. iHolman et ai]|2003[ ). A more direct con- 
nection was established between the observations and 
the ac celeration process first by lHamilton fc Petrosian! 
(1992), fitting to high spectral resolutio n but narrow 
band obser v ations (|Lin fc Schwartzlll987f ). and later by 
IPark et all (11997|). fitting to broad band observ ations 
(e.g. lMarschhauser et ai]|1994HDingus et al.lll994f ). This 
was done in the framework of stochastic acceleration 
(SA) by plasma waves or turbulence. 

However, it is preferable to obtain the X-ray radiating 
electron spectrum nonparame trically by some inv ersion 
techni ques first att e mpted bv lJohns fc Lin! 1 19921) . Re- 
cently, lHinaiFin (|2003D and iKontar et all (|2004 ) ap- 
plied regularized inversion techniques to obtain the radi- 
ating electron flux spectra from the sp atially integrate d 
photon spectra observed by RHESSI (jLin et all [20021 ) . 
This is an important advance but it gives the spectrum 
of the effective radiating electrons summed over the whole 
flare loop, but not the spectrum of the accelerated elec- 
trons. This difference arises because hi gh spatial res- 
olution observations, first from Yohkoh (Masuda ct al. 



1991 iPetrosian et al.ll200l and now from RHESSI (e.g. 
Liu et al.l 120031 ). have shown that, essentially for all 
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flares, in addition to the emission fr om the loop foot- 
points (FPs) (e.g. Hoyng ct al. 198lJ), there is substan- 
tial HXR emission from a region near the loop top (LT) . 
Thus, the total radiating electron spectrum is a complex 
combination of the accelerated electrons at the LT and 
those present in the FPs after having been modified by 
transport effects. 

It is therefore clear that separate inversion of the LT 
and FP photon spectra to electron spectra would provide 
more direct in formation on the ac celeration mechanism. 
More recently, iPiana et al.l (|2007t ) have applied the reg- 
ularized inversio n technique to the R HESSI data, in the 
Fourier domain (jHurford et al.l 12002) to obtain electron 
flux spectral images. The goal of this letter is to demon- 
strate that with the resulting spatially resolved electron 
flux spectra at the LT and FPs one can begin to constrain 
the acceleration model parameters directly. 

In the next section we present a brief review of the 
relation between the derived electron flux images and the 
characteristics of the SA model and in §3 we apply this 
relation to a flare observed by RHESSI. A brief summary 
and our conclusion are presented in §4. 

2. ACCELERATION AND RADIATION 

The observations of distinct LT and FP HXR emis- 
sions, with little or no emission from the legs of the 
loop, point to the LT as the acceleration site and 
require enhanced scatterin g of electrons in the LT. 
IPetrosian k, Do naghv (1999) showed that the most likely 
scattering agent is turbulence which can also accelerate 
particles stochastically. In fact SA of the background 
thermal plasma has been t he leading mechanism for ac- 
celeration of electr o ns (e.g. [Hamilton fc Petrosianlll992j: 
Miller et all 119961: IPark et all 119971: IPetrosian fe Liu! 
20041: iGrigis fc Benzl 120061: iBvkov fe Fleishman! 120091) 
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and ions (e.g. iRamatvl [19791: iMason et all 119861: 
Mazur et al]|1995t iLiu et al.ll2004 120061: iPetrosian et al.l 



2009), and is the most developed model in terms of com- 
paring with observations. 

2.1. Particle Kinetic Equation 

In this model one assumes that turbulence is pro- 
duced at or near the LT region (with background elec- 
tron density tt.lt, volume V, and size L). In presence 
of a sufficiently high density of turbulence the scat- 
tering can result in a mean scattering length or time 
(Tscat) smaller than L or the crossing time (t ctoss = 
L/v) , leading to a nearly is otropic pitch angle distribu- 
tion ([Petrosian fc Liull20M ). The general Fokker-Planck 
equation for the density spectrum N(E) of the acceler- 
ated electrons, averaged over the turbulent acceleration 
region, simplifies to 
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where Dee(E) and A(E) are the diffusion rate and direct 
acceleration rate by turbulenc^E respectively, El is the 
electron energy loss rate, and Q{E) and N(E)/T esc (E) 
describe the rate of injection of (thermal) particles and 
escape of the accelerated particles from the acceleration 
region. For electrons of energies below ~1 MeV, which 
are of interest here, Coulomb collisions^ dominate the 
energy loss rate, 



El = E Coul = 47rrgm e c 3 n LT lnA//3, 



(2) 



where In A is the Coulomb logarithm taken to be 20 
for so lar flare conditions. Following IPetrosian fe Liul 
(120041 ) . we approximate the escape time as T esc (E) ~ 
Tcross (1 + T cross /T scat ), which smoothly connects the two 
limiting cases of T cross /r scat 3> 1 and <C 1. The 
mean scat tering time is related to the pitch angle diffu- 
sion rates (jDung fe Petrosianll994l : |Pryadko fe Petrosian! 
[1993) due to both Coulomb collisions (D%° nl and r*^) 



and turbulence (£)*,„ and r s 



Tscat (E) = - 
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Similarly we can define the scattering times r s 



Coul 



and 



r turb 



for each process alone. For Coulomb colli- 
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For turbulence, D 



turb 



like Dee, depends on the spectrum of turbulence and 
on the background plasma density, composition , tern 
perature, and 



ma gnetic field (see [Sc hlickeiser 
jn[ [L994; Prvad ko fe Petrosian! 



1989; 



Dung & P etrosian! Il994t IPrvadko fc Petrosian! 11997 
1998L 119991: IPetrosian fe Liul 120041 ) . Since these coeffi 
cients determine the spectrum of the accelerated elec- 
trons, one can then constrain some aspects of the acceler- 
ation mechanism if an accurate spectrum of the electrons 
can be derived from observations. 

2 For stochastic acceleration, A(E) = D-eeC(E) / E + dDEE/d-E, 

where Q(E) = (2 - 7" 2 )/(l + 7" 1 ), 7 = l + E/m c c 2 = - f3 2 

is the Lorentz factor, and v = c/3 is the electron velocity. 

3 At higher energies, synchrotron loss must be included in E^. 



2.2. LT and FP Spectra 

The accelerated electrons in the (LT) acceleration re- 
gion with a flux spectrum Fm(E) — vN{E) produce thin 
target bremsstrahlung emissivity (photons s _1 kcV -1 ) 



J LT (e) = nurV / F LT (E)a(e, E)dE, 



(4) 



where o~( e, E) is the angle-ave raged bremsstrahlung cross 
section (jKoch fc MotJ[l959h . The escaping electrons 
with flux Fq(E) = N(E)L/T CSC produce thick target 
bre msstrahlung emi s sivity (coming m ostly from the FPs) 
(sec IPetrosianlll973l: IPark et alJfl997h . 

/OO 
F FP (E)o-{e 7 E)dE, (5) 

where n is the density and Fep is the effective radiating 
electron flux spectrum at the FPs, 



E L (n)J E T„(E') 

Since El oc n, the FP photon spectrum is independent 
of density. In what follows we evaluate equations ([5]) and 
^ using the LT density tilt- 

2.3. Acceleration Model Parameters 

Regularized inversion of RHESSI count vi sibilities 
gives the electron visibilities (jPiana et al.l 120071 ). which 
can then be used to construct images of electron flux 
(multiplied by column density). From these images, we 
extract the spatially resolved spectra, Flt(E) at the LT 
and Fep(E) at the FPs. Thus we can obtain the acceler- 
ated electron spectrum N(E) at the thin target LT. Also 
from differentiation of equation ([6]) we derive the escape 
time as T esc = — N(E)/-^(F F pE L /v), and by converting 
the denominator to a logarithm derivative we get 



Tesc (E) 



tl(E){F lt /F fp ) 
S FP (£)+2/( 7 + 7 2 ) 



t l (£)£(£), (7) 



where the FP index S FP (E) = - d , 1/(7 + 7 2 ) = 

~ d |iinE^ ' an< ^ t l(E) = E/El is the Coulomb loss time 
at the LT. The function £(E) is an observable quantity 
representing the ratio T csc /tl- In the above derivation, 
we have used the relativistic form of electron velocity 
v(E). 

Given T esc (E), from its relation to T cross and T sca t, we 
obtain the mean scattering time as r sca t — T 2 ross /(T csc — 
Tcross), which is valid for T C sc > T cross . Disentangle- 
ment of Ts" at b from T SC at is complicated (eq. [3]) at ener- 
gies when turbulence and Coulomb collisions contribute 
equally to r sca t- However, if turbulence dominates the 
pitch angle diffusion, then to the first order we can write 
T scat (1 + T SC at AS'), and obtain some average 



,_turb 



value of D 



turb 



Furthermore, given N(E) we can in 



principle determine the other Fokker-Planck coefficients, 
namely A(E) and Dee (see eq. PQ). Therefore we can 
reach a consistent picture of the acceleration process due 
to turbulence and begin to make inroads into the spec- 
trum and the nature of turbulence itself. 
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Figure 1. Electron flux spectral images (with 8 keV bin width above 34 keV and 2 keV bin width at lower energies) up to 250 keV 
in the 2003 November 03 flare during the n onthermal peak as reconstructed from two sets of the regularized electron visibilities by the 
MEM3IJIT algorithm (Schmahl ct al. 2007). The images show one LT and two FP sources above 34 keV and a loop structure at lower 
energies. Three circles are used to extract the LT and FP electron flux spectra above 34 keV (see Figure 
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Figure 2. Top: Electron power spectra MutE 2 F(E) for the LT (square), the two FPs summed (diamond), and all three sources (LT 
+ FPs, cross) in the 2003 November 03 flare. The LT spectrum can be fitted by a power-law, and the summed FP and total spectra by 
a broken power-law. Also note that the southern FP spectrum (downward triangular) is flatter than the northern FP spectrum (upward 
triangular), mostly above ~90 keV by ~0.3 powers of energy, consistent with their asymmetric locations with respect to the LT. Bottom: 
Escape time (filled circle) and turbulence scattering time (filled triangular) in the (LT) acceleration region. The escape time can be well 
fitted by either a power-law or a broken power-law (dash dot) increasing with energy, and the turbulence scattering time by a power-law 
rapidly decreasing with energy. Also shown are the crossing, Coulomb scattering, and the mean scattering (open triangular) times. The 
reduced chi-squares for all the fittings are below or around 1. 



3. APPLICATION: THE 2003 NOVEMBER 03 FLARE 

As a first demonstration, we apply our new procedure 
to the 2003 November 03 solar flare (X3.9 class) during 
the nonthermal peak, in which we find a hard LT source 
(extending above 100 keV in HXR) distinct from the 
thermal loop in addition to two FP sourcefl In Figure 
Q]we show the electron flux images up to 250 keV, which 
also show a loop at low energies and one LT and two FPs 
at higher energies. In Figure [5] top panel we show the 
electron spectra NurE 2 F(E), where Alt = nu^L is the 
LT column density. The LT flux spectrum can be fitted 
by a power-law with an index <5lt = 3.0. The summed 
FP flux spectrum can be better fitted by a broken power- 
law with the indexes <5i = 2.1 and 82 = 2.8 below and 
above the break energy E\> — 91 ± 3 keV. It is clear that 
the total radiating electron spectrum differs significantly 
from the (LT) accelerated electron spectrum. 

Given the above LT and FP electron flux spectra we 
derive the energy dependence of the escape time (eq. [7]). 
The LT density can be estimated as tt.lt — -JEM./L 3 ~ 
0.5 x 10 11 cm -3 , where the LT size L ~ 10 9 cm is ob- 
tained from the LT angular size, and the emission mea- 

4 Q. Chen & V. Petrosian (2010a, in preparation) present HXR 
observations of this flare and argue that the high energy LT source 
should not be an artifact of the pulse pileup effect. 



sure EM ~ 0.2 x 10 cm~ J is obtained from spectral 
fitting of the LT thermal emission. As in Figure [2] bot- 
tom panel, the escape time increases slowly with energy 
and can be fitted by either a power-law, 



T CSC (E) =0.3 s 



E 



100 keV 



k = 0.83 ±0.10, (8) 



or a broken power-law with a break at = 118 ± 37 
keV, and the indexes m = 0.62 ± 0.23 and K2 — 1-09 ± 
0.25. The fact that the escape time should be longer than 
the crossing time yields an upper limit on Alt > which is 
satisfied by the above LT density and size. 

We then calculate the mean scattering time in the LT 
region. Except at the lowest energy, the Coulomb contri- 
bution is small so that the scattering time thus calculated 
can be attributed to turbulence. The scattering time 
due to turbulence (see §2.3) can be fitted by a power-law 
above ^40 keV, 



turb 



0.016 s 



E 



100 keV 



A = 1.90 ±0.14. (9) 



4. SUMMARY AND DISCUSSION 

In this paper we describe a new method to directly 
obtain the model parameters for stochastic acceleration 
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of particles by turbulence in solar flares from regular- 
ized inversion of th e high resolution RHESSI HXR data 
(|Piana et al.l 120071) . We have argued that particle accel- 
eration takes place at or near the LT region. The accel- 
erated electrons produce thin target emission at the LT 
and then escape downward to the dense FP region un- 
dergoing Coulomb collisions and producing thick target 
emission. In this model the LT and FP electron spectra 
are connected by the escape process from the LT region 
(eq. [6]), thus allowing us to determine the energy de- 
pendence of the escape time. Our method has the ad- 
vantage that one can now constrain the model parame- 
ters uniquely rather than just satisfying the consistency 
between the model and the data as commonly done by 
forward fitting routines. This method can be applied to 
flares with simultaneous HXR emission from the LT and 
FP sources. 

We have applied our method to the 2003 November 03 
flare, in which we can obtain the electron flux images for 
both the LT and FPs up to 250 keV. The LT acceler- 
ated electron flux spectrum can be fitted by a power-law 
and the effective radiating flux spectrum at the FPs is 
better fitted by a broken power-law. From these spectra 
we derive the energy variation of the escape time and 
the scattering time. As seen in Figure^ the turbulence 
scattering time is relatively short and decreases with en- 
ergy. A short scattering time may arise from a high en- 
ergy density of turbulence (£turt>), with the exact rela- 
tionship depending also on the magnetic field (B), and 
the spectral index (q) and minimum wave number (fc m ; n ) 
of turbulence. A high level of turbulence also implies 
efficient acceleration which generally means a flat spec- 
trum for the accelerated electrons, which is the case for 
the current flare. The energy dependences of t*,^ and 
Dee are also a function of these characteristics of turbu- 
lence; at high energies they are deter mined primarily by 
the spectral index of turbulence ( see iDung fc Petrosian 
1991 IPrvadko fc Petrosianl H997L I1998L 119991: ILiu et al. 
2006). 

For the usually assumed Kolmogorov (q = 5/3) or 
Iroshnikov-Kraichnan (q — 3/2) turbulence spectra, one 
expects the scattering time to increase with energy as 
E 2 ~ q , which translates into an escape time varying 
roughly as T csc cx \j\J~E at high (but non-relativistic) 
energies. The energy dependences of T esc and r*"^ 13 ob- 
tained here require a steeper turbulence spectrum (q > 3) 
at high wave numbers. Such a steep spectrum can 
be present beyon d the inertial ran ge where damping 
is important (e.g. Jia ng et al.l [200 90 . The electron en- 
ergies and the wave-particle resonance condition deter- 
mine the wave vector of the accelerating plasma waves. 
This relation depen ds primarily on the pl asma parame- 
ter a oc \fnjB (e.g. iPetrosian fc Liu! 120041 ) . Thus, given 
the magnetic field and plasma density we can determine 
the wave vectors for transition from the inertial to the 
damping ranges of turbulence. 

It should, however, be emphasized that the results ob- 
tained here may not be representative of typical flares. 
More commonly flares have much softer LT emission, 
which would give an escape time decreasing (and scatter- 
ing time increasing) with energy, consistent with a low 
level and a flat spectrum of turbulence. 

The exact relation between the derived quantities 
(N(E), T osc , and t^?) and the turbulence character- 



istics (<?turb,-B,q, etc.) is complicated and depends on 
the angle of propagation of the plasma waves with re- 
spect to magnetic field and other plasma conditions. In 
future, we will apply these procedures to more flares (Q. 
Chen & V. Petrosian, 2010b, in preparation) and deal 
with these relations explicitly. 
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